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Abstract 

In the preceding paper, we derived Coulomb-gas and sine-Gordon Hamilto- 
nians to describe the Kosterlitz-Thouless transition on a fluctuating surface. 
These Hamiltonians contain couplings to Gaussian curvature not found in a 
rigid flat surface. In this paper, we derive renormalization-group recursion re- 
lations for the sine-Gordon model using field-theoretic techniques developed 

to study flat space problems. 
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I. THE SINE-GORDON HAMILTONIAN 



The Hamiltonian and associated partition function describing p-atic order on a fluctuat- 
ing surface were derived and analyzed in the Coulomb gas model in the previous paper [|l|. 
The renormalization-group (RG) recursion relations for K, k, and y were also derived from 
the Coulomb gas model. In addition, we described there how the two-dimensional Coulomb 
gas model can be transformed into a sine-Gordon field theory. 

Here we describe to what extent these results can be verified by a well-controlled renor- 
malization procedure based on the sine-Gordon field theory. The two-dimensional Coulomb 
gas model can be converted into the sine-Gordon Hamiltonian via a Hubbard- Stratonovich 
transformation. The advantage of this transformation is that it makes available to us stan- 
dard field theory diagrammatics and renormalization procedures |§||. This opens up a 
systematic way of obtaining results for the RG recursion relations. 

The p-atic membrane partition function in the sine-Gordon field theory is written as 

Z = J V^VRe-^-^-^^J^^, (1.1) 
where (3 is the inverse temperature, 

= / d 2 u^H 2 , (1.2) 

and 

^ = WK 7 rf2u vW*^aW- ^ j d 2 U^ COS0 (1.3) 

with = \K™ the mean curvature and S = det K* the Gaussian curvature. In the Monge 
gauge, the metric tensor g a p is written as 

f l + (d x h) 2 d x hd y h ^ 
. d x hd y h 1 + (d y h) 2 



g a /3 — d a Tl ■ dpR 
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and the curvature tensor K a p is 
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where (Vh) 2 = (d x h) 2 + (d y h) 2 . To lowest order in h, 



y d y d x h dydyh J 
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H = g a ^K aP = d x d x h + d y d y h = V' z h, 
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and 



S = det g af3 K aP = i(V 2 W 2 /i - didjhdidjh). 



(1.7) 

In performing a perturbation calculation with Hamiltonian (|1.2p and ( |1.3|) in two dimen- 
sions, one faces infrared as well as ultraviolet divergences. For the infrared regularization, 
we introduce a tension term 



[3H a = P<t[ d 2 Uy /g = P(t[ d 2 x (l + \(Vh)' 
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into the bending Hamiltonian ( |1.2| ) and a mass term 
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into the sine-Gordon Hamiltonian (|1.3j ). Our infrared treatment is similar to the infrared 
regularization in Ref. |J. After perturbation theory is summed, we shall see that our renor- 
malization group recursion relations have a well-defined limit when surface tension and mass 
are set to zero. Thus we will carry out our calculations in the presence of the tension a and 
the mass m, and set a —>■ and m — > at the end. 

The ultraviolet regularization is introduced by a short-distance cut-off, a, in the free 
propagators in coordinate space: 

47r 2 /3fr f d 2 q e iq - y 
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and 
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The Fourier transforms of these functions are 
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where i^oO^) an d ifi(x) are the conventional Bessel functions [Q and c = e 27 /4 with 7 
the Euler constant. We use this ultraviolet regularization rather than the sharp cut-off in 
momentum space since our calculations can be easily evaluated in coordinate space [§,[§,§1 • 

II. RG RECURSION RELATIONS 

In order to establish the RG recursion relations for the bending rigidity, the hexatic 
rigidity, and the fugacity of the fluctuating hexatic membranes, we study the renormaliza- 
tions of the two-point vertex functions T^(q) and T^l(q) for the total Hamiltonian in fll.l|). 
We have two independent sets of expansion parameters which are small. One consists of 
the single parameter, (/3k) _1 ; the other consists of the expansion parameters of the 2D XY 
model: the deviation of the hexatic rigidity from the fixed point value, 5 = (nf3K/2p 2 ) — 1, 
and the fugacity of the disclinations, y. We will carry out calculations to lowest order in 
temperature [(/?/t) -1 ] and combined second order in S and y. Thus we divide the renormal- 
ization scheme into two separate parts; first we obtain the vertex functions for y = 



up to first order in temperature, then we turn on the effect of the disclinations (y 0) and 
calculate the vertex functions r^ 2 ^ up to combined second order in S and y. 

~ (2) 

Starting with y — 0, we calculate the two-point vertex function for the sine-Gordon 

~ (2) 

field (j) to lowest order in temperature. The only diagram contributing to in this limit 
is shown in Fig. |IJ It arises from the coupling of to h. When a = 0, the resulting vertex 
can be written as 

f>($r \ ( P 2 p 2 3 1 \ 2 , P 2 2 (oU 

Thus, the effect of height fluctuation is to shift the hexatic rigidity from K to K, where to 
leading order in temperature 

{(3K)- 1 = my 1 + tJ-(^)- 2 - (2-2) 

OZ7T 

(2) 

Now taking into account the effect of the disclinations, we calculate up to second 

(2) 

order in o and y. The first step in the computation of rVi is the summation of tadpoles. By 
expanding the cos0 interaction of the sine-Gordon Hamiltonian in powers of 0, one gener- 
ates diagrams composed of vertices of all even orders and bare propagator G^^q, a). In Ref. 
0, Coleman noted that only diagrams containing tadpoles are ultraviolet-divergent loga- 
rithmically and these divergences can be removed by renormalizing the fugacity y. Tadpoles 
diagrams are shown in Fig. |2|. The effect of tadpoles is to renormalize y. Any diagram can 
be described as a diagram without tadpoles plus tadpoles adjointed. The effect of adding 
an arbitrary number of tadpoles to a bare vertex is to multiply each vertex by 

e -M*(*=°.«0, (2.3) 



thereby renormalizing y : 
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where the last equation defines A. 



Next we introduce some diagrammatic notation. The sum of even numbers of intermedi- 
ate propargators is [coshG^(x) — 1], and the sum of odd numbers of intermediate propar- 
gators is [sinhG^(x) — G^Jx)]. The diagrams corresponding to these sums are shown in 
Fig. |3| and Fig. |j. Using this notation, all diagrams contributing to up to second order 
in S and y are shown in Fig. [5]. Their sum can be written as 

- {yAf j d 2 x [e^[sinhGy*) ~ G%(x)] - [cosh G^ar) - 1]] . (2.5) 
To second order in S and y, the first order term in y contributes a factor of In a : 

yA = ycm 2 (cm 2 a 2 )~2p T = ycm 2 {l + 5\n(cm 2 a 2 )), (2-6) 



where 5 = (ir(3K/2p 2 ) — 1. The divergent part of the second order term in y has no q 2 - 
independent part. When we expand e tq ' x in powers of q, the odd terms in q vanish when 
integrated over angles. To second order in 5 and y, the coefficient of q 2 in is 

(yCTYl 2 } 2 f (yCTTl 2 } 2 r ml TTX^dx 

J d 2 xx 2 [smhG (x) - G {x)\ = ^ [cm2(a;2 + q2)]2+2 ^ + finite terms. 
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The integral contributes a factor of In a : 

pm~ 1 X^dx pm~ 2 zdz 

2(cm2) Vo [cm 2 {x 2 + a^ =L (7T^ = " ln(CmV) - (2 ' 8) 
The leading divergent contribution to is thus 
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rg(g) = -^^{(f + ™ 2 ) + %cm 2 [l + 5\n{cm 2 a 2 )} 

-^(2y) 2 q 2 ln(cm 2 a 2 ). (2.9) 

To establish the RG recursion relations for the hexatic rigidity (fiK)" 1 and the fugacity 
y, we define renormalized parameters (/3i^ r ) _1 and y r via 

{0KY 1 = Z K {f3Kr)- 1 , y = Z y y r . (2.10) 



All infinities can be absorbed by two independent renormalization constants. Zk and Z y 
can be chosen as functions of a to ensure that the renormalized two-point vertex, 



r i4, 2 j(9,2/rA,/i) = T {2) (q,y,5,a), 
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is finite for some arbitrary length scale fi, order by order in a double expansion in y r and 
5 r = (nf3K r /2p 2 ) — 1, in the limit a — > 0. Thus we expand Zk and Z y in a double series 
in y r and S r . These expansions are substituted into Eq. ( j2.10| ), which in turn are inserted 
in Eq. fl2.9|) . Then, terms of are rearranged by orders in the double expansion. The 
coefficient in Zk and Z y are determinied by the requirement that q 2 and m 2 have finite 
coefficients as follows : 
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The coefficients Zk and Z y are chosen so that Tp^l is finite in the limit a — > 0: 



, , 2 2ir 3 (3K r , 2 2 
Z^(/i) = 1 + ?/ r — ln/i a 
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where is an arbitrary mass scale. From the renormalization constants derived above, we 
obtain the RG recursion relations 
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and 
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where the subscript b means fixed bare parameters. In terms of length scale /x 1 = ae l , 
these equations are expressed as 



dl p 2 



{0K)- 1 = ~^V\ (2-17) 
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To complete the RG recursion relations for our Hamiltonian, we study the renormal- 
ization of r^-(g). We introduce a surface tension term [Eq. flOp ] to regularize infrared 
divergences. Strictly speaking, we should renormalize a HH as well as k. However, we will 
set the renormalized tension equal to zero in the end. We will, therefore, consider here only 
the q A contributions to T^, which will determine the renormalization of k. All diagrams 

(2) 

contributing to T hl i up to first order in temperature and second order in 5 and y are shown 
in Fig. ^| Their sum is written as 

rS(ff) = + ^ 4 M«m 2 ) - A^4 ln(cmV) 



3 16it 2 ({3K) 
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Defining renormalized parameter n r by 



3 7T 3 (^) 2 (22/) 2 g 4 (ln(cm 2 a 2 ))2 _ (219) 



k = Z K K r , (2.20) 



and using the renormalization constants Z# and Z y obtained in Eq. ( |2.13| ) and Eq. ( p,14[ ) 



(2) 

we find is finite in the limit a — > if we choose 



Z K (// = 1 - — — ln/i a +— — — — ln/x a 
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when expressed in terms of I = ln(yua) 1 



(2.23) 



Here we find y 2 term in Eq. (|2.21 ) is cancelled by the term proportional to y 2 from differ- 
entiating f3K r by fj, (See Eq. ( |2.16 )). Eq. Q2.23| ) does not have the term proportional to 



y ■ 

Thus the complete RG equations to first order in temperature and second order in 5 and 
y for our Hamiltonian are 



dl 



v 
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These equations and crikled-to-crumpled phase transition are analyzed in the preceding 
paper. 



III. ASYMMETRY IN THE DISCLINATION FUGACITIES 

On rigid flat membranes, a single positive disclination and a single negative disclina- 
tion have the same energy proportional to ln(R/a) where R is the linear dimension of the 
membrane and a is the size of the core region of the disclination. However, the plus-minus 
disclination energy symmetry present in rigid membranes is broken in deformable membranes 
|10| |. For deformable membranes, the topological charges q — ±1/6 of the disclinations can 
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partially be cancelled by the curvature charges due to the Gaussian curvature, when mem- 
branes buckle out to form a cone or a saddle surface. Membranes with a single disclination 
undergo a mechanical buckling transition from a flat configuration to a cone configuration for 
a positive disclination or a saddle configuration for a negative disclination. The Coulombic 
interaction between disclinations can be screened by Gaussian curvatures, and the disclina- 
tion energies are reduced in deformable membranes due to buckling into the third dimension. 
Furthermore, the energy of the cone is lower than that of the saddle, though both energies 



are proportional to \n(R/a) [pLl| . The system creates and destroys disclinations to minimize 
its free energy, and, in the hexatic phase, it will always choose configurations in which there 
is no overall \n(R/a) term in the energy. In flat membranes, this term can be eliminated 
by the condition of charge neutrality. In membranes whose shape can change, this term can 
also be eliminated by an average flat configuration and the condition of charge neutrality. 
We believe that this is the most likely scenario, however, it is not impossible that there 
are curved configurations with an overall unbalanced charge that manage to eliminate the 
\n(R/a) energy. Here we will consider only charge neutral configurations. Local buckling 
may lead to different fugacities y + and for positive and negative disclinations. However, 
because of the constraint of charge neutrality, physical observables such as the free energy 
will depend only on the product y+y~. 

We can generalize the our charge neutral renormalization calculation to take into account 
an asymmetry in positive and negative disclination fugacities. Using the Hamiltonian in Ref. 
IHII with different fugacities for positive and negative disclinations y + and respectively, 
we modify the sine-Gordon model as follows: 

Z = J V^VRe'^-^- 1 ^ I rf2 «v^ (3.27) 

where 

W, = ^(3k J d 2 u^gH 2 , (3.28) 

and 
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(3.29) 



Following the same renormalization scheme as above, we find the recursion relations for the 
bending rigidity, the hexatic rigidity, and the fugacities for positive and negative disclina- 
tions: 

d 47T 3 

j^kV = -rv + y-, (3.30) 



ir(3K 



*=|2-^fc (3.31) 
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Defining y = y/y+y~ and z = yy + /y_, we find 



mr 1 = Ary\ (3.33) 



d ( 7T(3K 



y=[2-^y, (3.34) 



Jf = °- (3 - 35) 
Equations for y and K are identical to symmetric case. Equation for z implies the ratio 
y+/y~ does not change under RG. Thus in the ordered phase, y + and y_ scale to zero and 
the KT transition occurs at the same temperature as in symmetric case. 

For T > T KT , both y + and y_ grow. If there is an initial asymmetry between y + 
and y_, it will persist, and positive disclinations will be favored over negative ones. This 
presumably favors positive Gaussian curvature. A complete discussion of the fluid phase 
would then require a more complete treatment of Gaussian curvature. The flat phase and 
the Kosterlitz-Thouless transition are, as we have just seen, not affected by an asymmetry 
in positive and negative disclination fugacities. 
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FIG. 1. Diagram contributing to the shift of the stiffness. 






FIG. 2. The circular dot represents the cos interaction and the square dot represents the sum 
of tadpole diagrams. 



FIG. 3. The sum of even numbers of intermediate propagators is represented by the vertical 
line of circles. 
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FIG. 4. The sum of odd numbers of intermediate propagators is represented by the sawtooth 



line. 
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FIG. 5. Diagrams contributing to T\i up to second order in y. The square dot represents 
the sum of tadpole diagrams in Fig. ||, the verticle line of circles the sum of even numbers of 
propagators shown in Fig. ||] and the sawtooth line the sum of odd numbers of propagators shown 
in Fig. U 






FIG. 6. Diagrams contributing to T\J up to one loop order and second order in y. The dotted 
line represents the height fluctuation field h. The square dot represents the sum of tadpole diagrams 
in Fig. the verticle line of circles the sum of even numbers of propagators shown in Fig. [|| and 
the sawtooth line the sum of odd numbers of propagators shown in Fig. 0. 
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